Photoemission spectra of a two-dimensional S=l/2 quantum antiferromagnet in 

magnetic fields: a theoretical study 
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We calculate the angular resolved photoemission spectra (ARPES) of a spin-1/2 Heisenberg an- 
tiferromagnet as a function of magnetic fields using both the exact diagonalization method and 
self-consistent Born approximation. Below the saturation field Be, strong scattering between spin 
waves and a hole, created by photoemission of an electron, significantly narrows the quasiparticle 
band that is characterized by the lowering of the quasiparticle energy at (tt, tt) with increasing field. 
Accordingly, in ARPES the quasiparticle peak gets shaper near (tt, tt) and broader elsewhere. Fur- 
thermore, we observe that an anomalous extended van Hove region (EVHR) around (7r, 7r) appears 
in a half saturation field, while EVHRs around (tt, 0) and (0, tt) in zero field gradually disappear 
with increasing field. 

PACS numbers: 75.50.Ee, 75.40.Gb, 79.60.-i 



The effect of a strong magnetic field on quantum spin 
systems has played an integral role in the understand- 
ing of magnetism and quantum phase transitions. There 
is renewed interest in this topic due to the synthesis of 
a new family of low-dimensional quantum antiferromag- 
nets with small exchange constants Q, 0, IH Q. For 
example, the theoretical prediction of field-induced in- 
commensurate soft modes in one-dimensional (ID) spin- 
1/2 quantum antiferromagnet (AFM) was confirmed by 
neutron scattering experiments on copper benzoate 
while anomalous spin excitation spectrum of a 2D AFM 
in strong fields has been reported [5J, l| ■ Other impor- 
tant observations include superconductivity in a layered 
organic AFM at very high magnetic fields 0, @ • In this 
paper, we present a theoretical study of the magnetic 
field dependence of the angular resolved photoemission 
spectra (ARPES) in a 2D AFM. 

Zero-field ARPES experiments in a 2D AFM are im- 
portant to the research of many cuprate superconduc- 
tors whose undoped parent compounds are nearly square 
lattice spin-1/2 antiferromagnetic insulators §. It was 
revealed that the presence of strong scattering processes 
between spin waves and a hole, created by emission of an 
electron, gives rise to an extended van Hove singularity 
near the Fermi surface. This novel feature is crucial in 
understanding many other anomalous physical properties 
of the cuprate superconductors and could be explained 
in the spin-polaron picture [ToL Hl| . On the other hand, 
above the critical field Be, the system becomes a satu- 
rated ferromagnet and the ARPES are known to be free- 
particle like. Therefore, it would be interesting to study 
the evolution of the ARPES of a 2D AFM from one limit 
to the other, in particular, the change of the extended van 
Hove regions. To our knowledge, this problem remains 
unanswered, primarily because the exchange constants in 
typical planar cuprates are J ~ 1500 K, hence, only near 
zero-field studies were possible. Nevertheless, a newly 
fabricated family of spin-1/2 square lattice antiferromag- 



nets, (5CAP) 2 CuA:4 and (5MAP) 2 CuX 4 with X=C1 or 
Br 0, 0, H, was found to have small J ~ 0.57 - 8.5 K 
with Be ~ 2 — 24 T, thus provides a good testing ground 
for our theoretical analysis. Alternatively, the ARPES 
studies could be carried out in certain pseudospin sys- 
tems where the effective magnetic field B e g is compara- 
ble with the effective exchange constant J e g; for instance, 
in a mixed valence system with electronic ferroelectric- 
ity 0, B e e is the d- and /-level energy difference 
that could be as large as J c g and adjustable by alloying 
or applying pressure. In the present work, we observe 
a magnetic-field-induced anomalous flat region around 
(tt, tt), where the bottom of the quasiparticle band lo- 
cates, in a magnetic field of \Bc, signifying an extended 
van Hove singularity near the Fermi surface from zero 
field up to \Bc- 

Our starting point is a 2D t-J model, a simple yet ef- 
fective Hamiltonian to model layered cuprate supercon- 
ductors, in a magnetic field along z direction, 



H 



(i,j> 



ia- c j<y + C j<x c i< 



where c icr = c i(T (l — n^) is the constrained fermion op- 
erator, and f ; = Y^vvAv?^™ with bein S the 
Pauli matrices is the spin operator. The SU(2) sym- 
metry is broken in the presence of an applied magnetic 
field, for the AFM will orient itself in such a way that 
the staggered direction, which is chosen to be the x di- 
rection, is perpendicular to the applied field. The result- 
ing magnetic phase is a canted state of two sublattices 
in which the spins tilt towards the z axis by the angle 
9 = arccos(£? / Be) with the saturation field Be = 4 J. 
Note that 8 = for B > Be- The Hamiltonian is invari- 
ant under rotation of spins around the z axis, hence the 
spin excitation spectrum should be a gapless Goldstonc 
mode. 
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The ARPES are defined as 

A(k,w) =J2 l(l^|ck,<r|0)| 2 5{u-Eo + Ei tV ), (2) 

where E and |0) are the ground-state energy and eigen- 
vector without any hole, and E lv and |1,^) are the en- 
ergy and the wave vector of the z^-th eigenstate with one 
hole created by photoemission of one electron. To calcu- 
late ^4(k, a;), we first employ the Lanczos exact diagonal- 
ization (ED) algorithm [14J with 100 iterations and an 
artificial broadening factor r\ = 0.05i on a 4 x 4 square 
cluster. Furthermore, in order to overcome the finite size 
effect in the ED calculations and gain more insights, we 
perform in the following analytical calculations of j4(k, lo) 
in the spin-polaron picture in the context of the self- 
consistent Born approximation. 

To simply the notation, it is convenient to perform a 
rotation of the spins in the A and B sublattices by 9 and 
— 9 about the y axis, respectively: 

~ ^ ; iOr ■ & i 

c iT = cos-d iT -e w 'sin-d u , 

Cil = e iQ r - sin ^d iT + cos ^d n , 

t? = S?cose + S?e iQr *sm6, 
t? = S? cos6 - Sfe iQ r * sin6, 

where tV — S\ and Q = (ir, ir). dia are constrained 

fermion operators and S i = Y2uv ^[uFiiv^iv are s P m °P~ 
erators in the new local coordinate system. This canon- 
ical transformation maps the canted spin configuration 
to a ferromagnetic configuration with all spins up and 
removes further necessity to distinguish between sublat- 
tices. The Hamiltonian thus has the form H = H t + Hj, 
where 

H t = -tcos6»^(4 T rf jT + 4x rf j|) 

-tsin0^e iQ - r '(4 T djj. -dld^ + H.c, (3) 
(id) 

Hj = jY,(S?S?cos29 + S?S?cos29 + S?S?) 

<iJ> 

+ J sin 29 (S? Sf - S? S? ) 

(id) 

+Bsm9j2s?e lQ r > ~ Bcos9j2s?. (4) 

i i 

Then, we treat quantum spin fluctuations within linear 
spin wave theory 0, 0] , namely, we regard the devia- 
tions of the spins measured from their equilibrium direc- 
tions are small, hence only up to the quadratic terms 
of the deviation operators are retained in Hj. Since 
the equilibrium directions are determined from vanish- 
ing of the linear terms, the second and third terms 



in Hj that are linear in spin deviations could be ne- 
glected. Note that there are spin-flipping hopping terms 
oc e l Q' ri {d\^d^ — in H t , which represent emitting 

and absorbing spin excitations, thus contribute to strong 
scattering processes between the hole and spin waves. 
Taking the ferromagnetic configuration as the vacuum 
state, we employ the slave-fermion formalism to cop e 
with the constraint of no doubly occupancy in H t |16| . 
Defining holon (spinless fermion) operators h\ so that 
= h\, dii = h\a i where eti = is the hard-core 
boson operator, we arrive at an effective spin-polaron 
Hamiltonian in the momentum space 

H - Y £k h i h k + Y w q aq« q 

k q 

+ ^(Mk q /jjA- q a q+Q + H.c), (5) 

k,q 

where a q 's are spin wave operators, a q = u q a q + t' q aL q , 
with dispersion w q = (A^ — 5 2 ) 1 / 2 . The transfor- 
mation coefficients are u q = [(A^/lu^ + l)/2] 1 ' 2 and 
u q = -sgn(B q )[(^ q /w q - l)/2] 1 / 2 . Here the shorthand 
notations are Aq = 2J(1 + cos 2 #7 q ), -B q = 2Jsin 2 #7 q , 
and 7 q = (cos q x a + cosq y a)/2 with the lattice constant 
a being the length unit in the following. Note that the 
spin- wave spectrum is indeed gapless (w q= Q = 0). The 
bare hole dispersion is e k = 4i cos 07k, and the hole-spin- 
wave coupling is 

At 

M kq = sin6i^=(7 k _ q w q+ Q + 7kV q+Q ). (6) 

v Jy 

Using the self-consistent Born approximation (SCBA) in 
which the spectral functions of one hole in a t-J-like 
model can be accurately calculated 0, 0, ITtI fisl llf| , 
that is, we first compute the hole Green's function 
G(k, uj) = [u> — e k — £(k,aj) + i0 + ] _1 self-consistently 
with the self-energy, 

E(k,c)=^M 2 q G(k-q,c-c q+Q ). (7) 
q 

Thus, the spectral functions of the hole quasiparticle 
(QP) are given by A(k,uj) = — ImG(k, lu)/tt, the spec- 
tral weights are Z(k) = [1 — <9£(k, uj)/dui\ ~= Ek , and the 
QP dispersion is E^ = + Re£(k, E^). 

To examine the finite size effect, we start the SCBA 
calculations with a 4 x 4 square lattice and compare the 
results with those obtained from using ED, then grad- 
ually increase the lattice size. Finally, we find that the 
SCBA results obtained from a 32 x 32 lattice are rather 
close (within 1%) to those from a 16 x 16 lattice. This 
implies that the SCBA results on the 16 x 16 lattice are 
not seriously affected by the boundary effect. Since in 
typical cuprates, J ~ 0.2 — 0.5t, we will adopt J = OAt 
in our calculations. The energy mesh is chosen from — 6t 
to 6t with interval O.Oli. 
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FIG. 1: The ARPES A(k,u) of the 4x4 antiferromagnet for 
J = OAt and B = \B C using (a) ED and (b) SCBA. (c)-(e) 
present the ARPES of the 16 x 16 antiferromagnet along (c) 
(0, 0)-(7T,7r), (d) (7r,7r)-(7r,0) and (e) (tt,0)-(0, 0) directions 
using SCBA. 
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FIG. 2: Same as in Fig.EJ except for B = f Be- 



Then we proceed to compare the ED and SCBA re- 
sults. From the ARPES A(k, lu) of the 4 x 4 AFM in 
Figs. IHa)-(b) for B = \B C and in Figs. E{a)-(b) for 
B = \Bc-, we notice that the line shapes obtained from 
using ED [Figs.^a) and^a)] agree well with those ob- 
tained from SCBA [Figs. □» and l^b)]. As for the 
results for B = 0, we refer readers to Ref. 0, Il7j . 
These numerics unambiguously demonstrate that the 
spin-polaron picture provides indeed a natural descrip- 
tion of the QP behavior. In addition, the ARPES of 
the 16 x 16 AFM obtained from using SCBA shown in 
Figs, dc)-(e) and Fig. Hlc)-(e) are in consistence with 
those obtained from the 4x4 AFM. The main feature 
of these ARPES is that as B increases, the spectrum at 
(-7T, 7r), the QP band bottom, gets more and more coher- 
ent, while away from the band bottom the spectra be- 
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FIG. 3: Quasiparticle dispersion for J = OAt and B = (a), 
\Bc (b), and \Bc (c). Left label: ED results (solid squares) 
on the 4x4 square lattice. Right label: SCBA on the 4x4 
(open circles) and 16 x 16 (solid lines) lattices. 



comes incoherent quickly, indicating that the QP states 
away from the minimum decay by emission of spin waves. 

Furthermore, we present in Fig. the calculated elec- 
tronic structure of the 2D Heisenberg AFM with B = 0, 
\Bq and |-Bc, respectively. Again, the SCBA results 
agree well with those obtained from ED. We notice three 
main features in these quasiparticle bands: First, as B 
increases, the band bottom evolves from (tt/2,tt/2) to 
(it, 7t). Second, the low-field band [Fig. Of a)] has flat re- 
gions around (ir, 0) and (0,7r), while the high-field band 
[Fig. Ofc)] resembles the free particle dispersion but its 
width is severely narrowed. In between, for B = \Bc 
[Fig- Bl^h)], a new anomalous feature, a flat region around 
(•7T, 7r), appears in the electronic structure, which leads to 
a strongly distorted density of states with a massive peak 
near the bottom of the QP band. Therefore, an extended 
van Hove singularity near the bottom of the QP band 
survives up to B — \Bq. Third, the QP bandwidth is 
strikingly narrowed, which can be seen more clearly in 
Tabled 

In Table[IJ we present several interesting physical quan- 
tities, including the QP bandwidth W, the spin-wave 
bandwidth W sw , the minimal QP energy E mm , and the 
QP spectral weights Z(k) at high symmetric points in 
the first Brillouin zone, as a function of the magnetic 
field. These data are obtained from using SCBA on the 
16 x 16 lattice. We find that W and W sw have similar 
values. Such band narrowing can be understood as fol- 
lows: since gapless spin excitations are easily stimulated 
by the incoherent motion of a hole, the combination of 
the hole and polarized spin-wave cloud constitutes the 
quasiparticle, spin-polaron. Hence, W scales with W sw . 
Another special feature revealed from Table |U is that in 
low fields the largest spectral weights locate at (tt/2, 7t/2) 
and (-7T, 0); increasing B would raise Z(tt, 7t), while it de- 
creases the QP spectral weights at other wave vectors. 
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TABLE I: -E min , the minimal QP energy, W, the QP band- 
width, W BW , the spin- wave bandwidth, and Z(k), the QP 
spectral weights at several high symmetric points in the first 
Brillouin zone as a function of B, the magnetic field, obtained 
on the 16 x 16 lattice using SCBA for J = OAt. 

B/Bc £min/t W/t W sw /t Z(0,0) gCLj[j Z{n,n) Z(tt, 0) 
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Heisenberg antiferromagnet. Using the self-consistent 
Born approximation in the spin-polaron picture, we ob- 
serve an anomalous electronic structure in which an ex- 
tended van Hove region around {%, 0) in zero field moves 
toward (n, tt) with increasing field strength until reaching 
a half saturation field. We also present the exact diago- 
nalization data that confirm our prediction, which could 
be tested by ARPES experiments on weakly interacting 
spin- 1/2 square lattice antiferromagnets sag. 

We are grateful to C. D. Batista and H.-Q. Lin for 
useful discussions. This work was supported by the Ne- 
braska EPSCoR-NSF Grant EPS-9720643 and the Office 
of Naval Research. 



This confirms the spectral features shown in Figs. ^ and 
121 We explain this feature as follows: in low fields the 
spins are almost antiferromagnetically aligned, therefore 
effective hole hopping in the same sublattice would be 
favored, leading to -E m in at (7r/2,7r/2) and considerably 
large values of Z(tt/2,tt/2) and Z(it,0). On the con- 
trary, increasing B will raise the bare hole dispersion 
e k oc cos 9 — B /Be and decrease the hole-spin- wave cou- 
pling Mk q oc sin#, thus facilitates the hole propagation 
to nearest neighboring sites and lowers the QP energy at 
k = (tt, 7r), whereas hole motion with other wave vectors 
is impeded due to scattering off spin waves. We notice 
from Table |I| that W and Z(k) experience abrupt change 
as B — > Be- In this limit, Eq. Q can be solved analyti- 
cally in perturbation theory, since Mkq — > as B — > Be- 
Hence, the QP dispersion is given by 



£k 



-E 



£k + 0(l- 



B_ 

B C ' 



(8) 



Above Be, Mk q = 0, thus — = 4t7k, that is, the 
hole moves freely in a saturated ferromagnet. 

It is worth mentioning that linear spin- wave theory was 
pointed out to be inaccurate in strong fields - spin waves 
in a quantum AFM in a strong field (> 0.765(7 ) were 
unstable with respect to the spontaneous two-magnon 
decays @. Nevertheless, we have demonstrated that 
the strong-field ARPES obtained from using the SCBA, 
where the spin excitations are treated in linear spin- wave 
theory, agree well with the ED results. Hence, it appears 
that linear spin-wave theory is appropriate to calculate 
the spectral functions of a hole in a quantum AFM. While 
it is desirable to follow the present study with high order 
spin-wave theory or ED for larger lattices, we anticipate 
that our basic conclusions remain valid. 

To summarize, we have studied the magnetic field de- 
pendence of the ARPES of a square lattice spin-1/2 
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